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Abstract. We present a model for the equilibrium of solid planetary cores embedded in a gaseous nebula. From this 
model we are able to extract an idealized roadmap of all hydrostatic states of the isothermal protoplanets. The 
complete classification of the isothermal protoplanetary equilibria should improve the understanding of the general 
problem of giant planet formation, within the framework of the nucleated instability hypothesis. We approximate 
the protoplanet as a spherically symmetric, isothermal, self-gravitating classical ideal gas envelope in equilibrium, 
around a rigid body of given mass and density, with the gaseous envelope required to fill the Hill-sphere. Starting 
only with a core of given mass and an envelope gas density at the core surface, the equilibria are calculated without 
prescribing the total protoplanetary mass or nebula density. In this way, a variety of hydrostatic core-envelope 
equilibria has been obtained. Two types of envelope equilibria can be distinguished: uniform equilibrium, were 
the density of the envelope gas drops approximately an order of magnitude as the radial distance increases to 
the outer boundary, and compact equilibrium, having a small but very dense gas layer wrapped around the core 
and very low, exponentially decreasing gas density further out. The effect of the envelope mass on the planetary 
gravitational potential further discriminates the models into the self-gravitating and the non-self gravitating ones. 
The static critical core masses of the protoplanets for the typical orbits of 1, 5.2, and 30 AU, around a parent star 
of 1 solar mass (Mq) are found to be 0.1524, 0.0948, and 0.0335 Earth masses (M^), respectively, for standard 
nebula conditions (Kusaka et al. I1970II . These values are much lower than currently admitted ones primarily 
because our model is isothermal and the envelope is in thermal equilibrium with the nebula. Our solutions show 
a wide range of possible envelopes. For a given core, multiple solutions (at least two) are found to fit into the 
same nebula. Some of those solutions posses equal envelope mass. This variety is a consequence of the envelope's 
self-gravity. We extend the concept of the static critical core mass to the local and global critical core mass. Above 
the global critical mass, only compact solutions exist. We conclude that the 'global static critical core mass' marks 
the meeting point of all four qualitatively different envelope regions. 

Key words, planetary systems: formation, protoplanetary disk; Planets and satellites: general; Solar system: 
general 



1. Introduction 

With the discovery of the extra-solar gas giants, the gen- 
eral problem of planet formation has considerably grown 
in complexity over the last decade. However, a global the- 
oretical overview of the properties of the giant planets, 
irrespective of the parent protoplanetary disc or the total 
mass of the giant planet, is still missing. 
In the nucleated instability hypothesis, envelopes of giant 
planets are thought to be formed as a consequence of ac- 
cretion of solid bodies forming their cores. To determine 
the envelope mass corresponding to a given core, static 
protoplanetary models have been constructed (e.g. Perri 
& Cameron [WT^ Mizuno [1^501 Stevenson . 1982,1 . 



If the envelope is modelled including detailed energy trans- 
fer and if the outer part of the envelope is radiative, and 
for standard assumptions about nebula conditions, it has 
been found that there is an upper limit for the masses of 
static envelopes and therefore for the total mass of a proto 
giant planet. This upper limit in core mass - the critical 
mass - was found to be insensitive to nebula conditions, 
but to depend weakly on dust opacities (Mizuno I1980|l 
and on the rate at which the core (solid body) is accreted 
(Stevenson Unil. 

Even the largest static critical masses are typically more 
than a factor of ten smaller than Jupiter's mass (Mizuno 
[TMil StevensonUnia Wuchterl lTMTbl Ikoma et aL lMmi . 
The nondependence of the critical mass on nebula condi- 
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tions disappears when the outermost parts of the proto- 
planetary envelopes become convective, which happens for 
nebula properties which are well within of proposed solar 
nebula conditions (Wuchterl I1993|) . Envelope masses of 
such protoplanets range between 6 and 48 Earth masses 
(Mq) but hydrostatic models alone are unable to re- 
produce a Jupiter-mass protoplanet. Therefore dynamical 
and/or quasi- hydrostatical effects should play an impor- 
tant role in the formation of gas giants. 
There is a number of incompletely studied processes (e.g. 
the formation, evolution, and stability of the protoplane- 
tary disks, the dust growth, the planetesimal formation, 
etc. ) that are relevant for the general problem of planet 
formation. Their complexity makes a piecewise approach 
necessary in studies of planet formation. An alternative 
approach is to study the final outcome, i.e. the possible 
and probable end-states of the process. In that context, 
we present an idealized road-map of all hydrostatic states, 
in order to provide insight when analyzing the complex 
behavior of hydrodynamic and quasi-hydrostatic models 
with detailed microphysics. In addition, this work aims to 
clarify the concept of the critical core mass necessary to 
permanently attract gas of the protoplanetary nebula to 
a terrestrial-planet-like heavy element core. 



2. Model 

2.1. Motivation 

Within nucleated instability theory, the formation of gi- 
ant planets includes many possible scenarios for pro- 
toplanetary cores and their respective envelopes. These 
range from small planetoids embedded in dilute proto- 
planetary nebulae to present-day-Jovian-like cores of sev- 
eral M0 squeezed by some Mbars of metallic Il2-IIe mix- 
tures (Guillot I1999|l . To date, many investigations have 
been made into the evolution of protoplanets, both hy- 
drostatically (e.g. Bodenheimer et al. I2UUUL Ikoma ODUll 
see Wuchterl et al. '2000| for review) and hydrodynami- 
cally (e.g. Wuchterl 1 199Tal I1991bl IT^Mjl . In these studies, 
'the evolution' of particular planets is followed, but not 
much is known about the evolution of all possible pro- 
toplanets. Therefore, it is somewhat difficult to bring the 
detailed solutions of previous investigations within a global 
perspective. 

We follow the thermodynamical approach that was used 
by Stahl et al. I)1995|l to investigate the coreless equilibria 
of constant-mass isothermal gas spheres, and the nature 
of the Jeans instability. We also expand on the work of 
Sasaki IT9891 who studied isothermal protoplanets in the 
minimum mass solar nebula (MMSN). In our model the 
total mass of the protoplanet and the density of nebula 
cloud in which the protoplanet is embedded are not pre- 
scribed. In leaving these as output variables, and starting 
only with the (heavy-element) core mass and the density 
of the envelope gas at the core's surface, we aim for a 
complete classification of all hydrostatic equilibria. This 
classification should contribute to clarifying whether mul- 



tiple planetary equilibria exist for given nebula conditions 
and how protoplanetary models relate to gas giants, both 
inside and outside of the solar system. 

2.2. Model Assumptions 

We approximate the protoplanet as a spherically symmet- 
ric, isothermal, self-gravitating classical ideal gas envelope 
in equilibrium around a core of given mass. This gaseous 
envelope is that required to fill the gravitational sphere of 
infiuence, approximated by the Hill-sphere: 



r-HiU 



a^Afpianot/3M^, 



(1) 



where a is the orbital distance from a parent star. With 
mean molecular weig ht of /X = 2.3 10-3 kgmol-\ pro- 
toplanetary envelopes, as well as the nebula, are roughly 
approximated by a hydrogen-helium mixture. The pro- 
toplanet 's heavy-element-core is represented by a rigid 
sphere of uniform density of Qcotc — 5500 kg m'^. 

The nebula temperature profile is taken according to 
Kusaka et al. flWOI, and Hayashi et al. l(Tn55jl . cf. TableHJ 
The nebula density structure is not a priori determined, 
but, for critical core mass determination, nebula densities 
agree with those from Kusaka et al. H197 0"l for a = 1 and 
30 AU, and from Hayashi for a = 5.2 AU, cf. Table 

12 It has been shown that the critical core mass values 
have only a weak dependence on the nebula density (cf. 
Sect. l3.9|l . therefore the choice of the nebula density is not 
critical. 

2.3. Model Equations 

The envelope is set in isothermal hydrostatic equilibrium, 
with spherical symmetry, and as such is described by: 



dM{r) 
dr 



= 47rr^p(r), 



the equation of hydrostatic equilibrium: 
dP{r) _ GM{r) 



dr 



and the equation of state for an ideal gas: 
P{r) = —g{r). 



(2) 



(3) 



(4) 



M(r) is defined as the total mass (core plus envelope) 
contained within the radius r: 



M(r) = Mcore + / 47rr'2£.(r') dr' 



(5) 



where r is the radial distance measured from the core cen- 
ter and g is the envelope gas density at radial distance 



Table 1. Symbols 
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Symbol^ 


Meaning 


Orb. param. (a, T) 


(1, 225) 


(5.2, 123) 


(30, 51.1) 


a [AU] 

G = 6.67259 10"" 
= 2.3 10"^ 


orbital distance 
gravitational constant 
mean molecular weight 


MSr/[Mffi] 


0.1524 
21 


0.0948 
96 


0.0335 
380 











,30 



Moore 
Menv 

Mtot 
M(r) 

Mq = 1.989 10' 
M® = 5.976 10^'' 

^corc 

rmii 

K = 8.31441 
fcorc = 5500 

Ocsg 
Qcnv 

T{a) 



predefined core mass 
envelope mass 
total mass 

total mass interior to radius r 

solar mass 

Earth mass 

core radius 

Hill sphere radius 

molar gas constant 

predefined core density 

envelope gas density at core surface 

envelope gas density 

nebula gas temperature 



t 



SI units used unless otherwise specified 



2.4. Boundary Conditions 

The total mass of the protoplanet is defined as: 



with 



Mr. 



Mr, 



M(rout) 



(6) 



A^(»-corc) =Mcore. (7) 

The inner and outer radial boundaries are: 



Mr, 



'in — ' core 



4 

3 ^QcOTC 



and Tout = fHiii- 



(8) 



An additional boundary condition at the core surface is: 

^?env (^coic ) Qcsg- (9) 

This model, together with the specified assumptions 
and boundary conditions, is sufficient to completely de- 
termine a single model- protoplanet. The total mass and 
nebula density at rniu (gas density at protoplanet 's outer 
boundary) are results of the calculation. 

2.5. Solution Procedure 

The total protoplanetary mass is obtained by integrating 
outward from rcore to r'Hiii(Mtot), starting with r^jjj = 
rHiii(Mcore) and iterating rHiii(Mcorc + Mmv)- 
Integration is performed from the core surface to the Hill 
radius, using the Maple 6 software (e.g. Garvan 2001), 
with the Fehlberg fourth-fifth order Runge-Kutta method. 



The critical core mass increases for smaller orbital distances 
because of (in order of importance): the higher gas temper- 
ature (cf. Sect. 13. 8"n and 15. 101 . the smaller Hill sphere (cf. 
SectE 



3, and the higher densities of the reference nebu- 
lae (taken from the minimum mass solar nebula models of 
Kusaka itTTTOl and Hayashi lITOSIl '). 



3. Results 

3.1. Solution Manifold 

In order to cover as many hydrostatic solutions as possible, 
the system of equations , © , and Q is solved for a wide 
range of parameters Mcorc and Pcsg- The set of all solutions 
for this range constitutes the solution manifold. Figure ^ 
shows the solution manifold for a protoplanet whose or- 
bital distance corresponds to the position of proto- Jupiter 
according to the Kyoto-model of solar system formation 
(Hayashi et al. 1985^) . The manifolds with orbital parame- 
ters (a, T) of proto-Neptune and proto-Earth have similar 
morphologies. It should be reiterated that the solution set 
contains all qualitatively different protoplanetary models 
at a particular orbital distance; not just for a particular 
nebula, but for any nebula, from a dense gravitationally- 
just-stable cloud to a near-vacuum space. 

3.2. Manifold Regions 

Several distinct regions exist in the parameter space of 
the solution manifold (Fig. 12), and they can be examined 
from two complementary perspectives. One way is to use 
the gas density at the core surface, Qcsg, as an independent 
variable (e.g. Fig. |2J), and the other is to use the nebula 
gas density, gout (e.g. Fig.^. While gout is more physically 
intuitive, Pcsg maps out region IV of Fig. |2] more clearly, 
and is more efficient in terms of representing the entire 
manifold. 

Figure [21 divides the solution manifold into four distinct 
regions, depending on whether the solution is compact or 
uniform and self-gravitating or not. Figures |2| and ^ point 
to the existence of the four possible regimes for a planet; 

1. 'mature telluric planet' (region I): envelope mass is a 
linear function of (Pout, and Qcsg 

2. 'mature giant planet ' (region II) : envelope mass weakly 
drops with gout (Monv oc £iout"°^)- 

A/cnv Of Qrsg is weaker than for the 'nebula' regime. 
'Nebula' densities (giout) are so low that they may well 
be considered a vacuum. 
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Fig. 1. Envelope mass solution manifold. Environmental parameters for this manifold are set to a=5.2 AU, and T—12S 
K. Each point on the surface gives the mass of the protoplanet's envelope for a given Afcorc and gas density at the core 
surface, Qcsg- Models with different initial parameters generally connect to different nebulae. Several different regions 
are easily discernible: I - flat slope with gradient of 1, for the region [-1,2] in logMcore and [-12,6] in log^pcsg; H - flat 
slope with gradient of 0.5, roughly encompassing [4-6,8] in loggcsg, and all logMcorc; HI - 'base of the island', [-8,-1] 
in logMcore and [-12,-6] in logfeg; IV - 'island', [-8,-1] in logMcore and [-6,4-8] in logfeg (cf. Fig.E)). 



3. 'ne^wZa '(region III): envelope mass is a linear function 3.3. Self-Gravity Effect 

of Pout and Qcsg 

4. 'protoplanet' (region IV): envelope mass is a non- 
trivial function of pout or Qcsf 



The key effect, which is responsible for the manifold mor- 
phology as observed in Fig. ^ can be described as self- 
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Fig. 2. Manifold regions: I - compact non-self-gravitating 
envelopes, II - compact self-gravitating envelopes, III - 
uniform non-self-gravitating envelopes, IV - uniform self- 
gravitating envelopes. The border of the region IV some- 
what depends on the choice of the surrounding nebula (cf. 
Fig ll2|l : we use here a value from the Hayashi (4985^) min- 
imum mass solar nebula model. 



gravity of the protoplanet's envelope. Keeping in mind the 
hydrostatics of the model, and the fact that the surround- 
ing nebula is not prescribed, one can see that self-gravity 
reduces the envelope mass for a given core surface pres- 
sure, i.e. the envelope mass would be larger if there were 
no self-gravitating effect (Fig. |3J|. 

The envelope's radial gas density profile is shaped 
through the interplay of inward gravitational force and 
outward gas pressure. If the envelope mass is small com- 
pared to the core mass, the gravitational force can be ap- 
proximated as arising from the core's gravitational poten- 
tial only. Once the envelope mass is comparable to (or 
greater than) the core mass, they both contribute to the 
gravitational potential, making its gradient steeper and, 
in effect, reducing the envelope mass. As a consequence, 
the self-gravitating envelope connects to a nebula different 
from the one which is in balance with the envelope in the 
absence of the self-gravitating effect. Further discussion of 
the role of self-gravity can be found in Sect. 13.51 

3.4. Two types of envelope equilibria 

The solution manifold (Fig.|2Jl contains two basic types of 
envelope equilibria (Fig.EJ: 

1. uniform, or quasi-homogenous envelope: the density of 
the envelope gas drops weakly with increasing radial 
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Fig. 3. Demonstration of the self-gravitating effect for 
sub- and super-critical cores: comparison of cuts through 
two manifolds, with- {M ~ M{r) in Eq. and without- 
(M = Mcore) the envelope's gravitating effect, each for 
two core masses. Cuts are for a = 30 AU and T = 51.1 
K. Circles and squares represent the envelope mass of the 
subcritical core, calculated for M = M{r) and M — Mcore 
in Eq. (j3Jl , respectively. White and black triangles have the 
same meaning but for the supercritical core. Labels with- 
out arrows correspond to manifold regions from Fig. 12 
while labels with arrows mark interfaces between regions. 
D corresponds to the 'divergent wall' which surrounds re- 
gion IV (cf. Fig^. Self-gravitating envelopes with M = 
M{r) in Eq. ^ have a larger envelope mass than the 
corresponding envelopes with M = Mcore in Eq. Q (cf. 
Fig. El. 



distance, keeping the mass distribution more or less 
uniform throughout the envelope; 9Monv/3rout > 
2. condensed, or quasi-compact envelope: typically small, 
but very dense gas layer is wrapped around the 
core, at larger radii the gas density is very low; 
9Mc„v/5rout » 

This is reminiscent of a similar equilibrium, found by Stahl 
I1995II. for constant mass coreless 'Van der Waals' 



et al. 

gas spheres. 

If an envelope's mass is much smaller than the core mass, 
the radial profile of the gas density is simply an exponen- 
tial function, well approximated by : 



P{r) = P,eM-^GM{r){^ - -)). 

j\l ^corc ^ 



(10) 



If (M(r) - Mcore) < Mc, 
barometric formula. 



then Eq. H10|l reduces to the 
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Fig. 4. Envelope mass as a function of the nebula density 
Qout- Labels are the same as in Fig.|21 Lines connect states 
with increasing gcsg- Note the strong dependence of gout 
on the envelope mass, and a non-trivial behavior of the 
McnviQout) for region IV (enlargement in Fig.|Sl). 

3.5. Differences: Isothermal Careless Gas Spheres vs. 
Isothermal Protoplanets 

The curl in Figs. 01 and El is reminiscent of a similar fea- 
ture found for the isothermal coreless ideal-gas spheres 
(e.g. Schwarzschild. [1^571 Sect. 13) represented in the U- 
V plane. It follows from the definition of U and V that: 

J J ^ r dM{r) ^ 47rr^g ^ g 

M{r) dr M{r) M{r)/(^TTr^) ^ ' 

y_ rdP^g GM{r) ^ 3 GM{r)/r 

P dr P r 2 |P/g ^ ' 

and from the fact that the mean density of the total object 
for our model is always the same, as implicitly defined 
through Eq. (jlj. 

Unlike a singular isothermal sphere with an infinite 
pressure at the center, our protoplanetary model has a 
solid core of miiform (and finite) density at its center. 
This will result in the departure from the potential of the 
coreless isothermal sphere: instead of having a g(r) oc 
structure, the envelope gas close to the core surface will 
obey the form of the barometric law (c.f. Eq. 
If the mean envelope density at lower stratifications is 
comparable to the core density, an 'effective' core will 
shorten the characteristic length-scale of the potential, 
making the exponential profile of the barometric-law-like 
profile even steeper. For the appropriate effective core, the 
outer stratifications will exactly match the outer stratifica- 
tions of the solution which has the gas density at the core 
surface much smaller than the core density (cf. Fig. 1141) . 
These profiles will connect to the same nebula density, 
but will have a slightly different envelope mass, because 




a&:H3"'9.I3)(H3"' l.CwHJ* t. lie IC* t.ailC* t,3x10'* 1 .'fxlO'* 1 .ftcIO'' I.SclO'* 

Fig. 5. Enlargement of the boxed region of Fig.01 isother- 
mal curl regularized with the finite-density core; '-1.25': 
black squares represent protoplanets with first subcritical 
Mcorc hne on the mesh of Fig.^and the arrow points at the 
black square with the highest Monv, DS: two protoplan- 
etary states with the largest envelope mass in the mani- 
fold, but with typically very different gcsg (cf. Sect. l3.9.'2)l : 
in and out curves are the consequence of the core. The 
smaller the core, the closer the in and out curves are. The 
figure corresponds to a V-U plane for the protoplanets (see 
Sect. 13.51 for further discussion). 



of the difference in the profile of the inner stratification. 
Therefore, the curl of Fig. O will have two branches: 'in' 
(the solution with a non-sclf-gravitating inner stratifica- 
tion) and 'out' (the solution with an effective core). 
The smaller the core mass, the sooner will the profile con- 
nect to the 'r~^' structure; i.e. the smaller the difference 
in the envelope mass between the pairs of solutions, the 
closer the 'in' and 'out' branches in Fig. will be. 

The fall-off of the gas density with increasing radius 
in the self-gravitating part of the envelope can be approx- 
imated by Qcnvir) ~ r~'^ (cf. Fig. self-gravitating pro- 
file) , as expected in the theory of stellar structure for a self- 
gravitating isothermal sphere of ideal gas (e.g. Shu "1 9921 
Sect. 18). Small deviations from are due to the finite 
amount of mass needed for the envelope to become self- 
gravitating, which produces a slight imbalance between 
the self-gravity and the amount of mass M(r). No similar 
effect is observed for coreless, isothermal gas spheres (cf. 
Stahl et al. iro^ . 

Depending on the fraction of the self-gravitating part of 
the envelope and of the core mass, this wavelike deviation 
can extend to the outer boundary, or can be attenuated 
deep within the envelope. 
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Fig. 6. Uniform, compact and self-gravitating profiles. 
The uniform self- gravitating profile resembles the non- 
self-gravitating one until the envelope mass becomes com- 
parable to the core mass. Then the density profile changes 
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Fig. 7. Pressure as a function of density, for T ~ 123 K. 
Black circles represent the ideal gas, squares are for the 
Carnahan-Starling EOS, and triangles are for the Saumon- 
Chabrier EOS. This figure also shows that a completely 
degenerate electron gas (stars) is not a good assumption 
for this {g, T) parameter range. 



3.6. Estimating the Applicability of the Ideal Gas 

We made two major assumptions while constructing our 
model - that the gas is ideal, and that the heat is in- 
stantaneously radiated away, i.e. the gas is isothermal. In 
Sect, im we examine the isothermal assumption, and we 
deal with the ideal gas in this section. 

In order to keep the protoplanet in an equilibrium with 
the surrounding nebula, we have set the envelope gas tem- 
perature equal to the nebula temperature for the appro- 
priate orbital distance. Therefore, we compare different 
equations of state at the envelope temperature. In addi- 
tion to an ideal gas, we take the Saumon-Chabrier-van 
Horn EOS, the Carnahan-Starhng l(TU^ EOS, as 

well as a completely degenerate electron gas. 

Figure shows that for the gas densities up to about 
40 kgm~'^, ideal gas, Saumon-Chabrier-van Horn, and 
Carnahan-Starling EOS agree to better than one percent. 
For higher densities the Saumon-Chabricr EOS shows 
additional non-ideal effects, while the Carnahan-Starling 
EOS exhibits a similar behavior for densities larger than 
200 kgm~'^. We can also see that the electron degeneracy 
does not contribute to the pressure at least till the point 
where the Saumon-Chabrier EOS departs from ideal-gas 
behavior. 

However, in general we see that the ideal gas is an ex- 
cellent approximation for our model for the better part 
of the envelope gas density range. Certainly, there are 
also models where densities are high enough for signifi- 
cant non-ideal effects, but typically for the protoplanets 
in our model those high density envelope regions are re- 
stricted to areas close to the core, while the rest of the 
envelope will be well approximated by an ideal gas. We 
can see in Fig.|Hlthat if we use e.g. the Carnahan-Starling 



EOS, the numerical details will be changed, but the qual- 
itative picture will remain the same. This is also true for 
the Saumon-Chabrier EOS, which is work in preparation 
by C. Broe.g. The ideal isothermal gas will not be a good 
approximation for the compact envelopes which are typ- 
ically associated with giant planets in the late stages of 
their evolution. Using our model, we can show that a pro- 
toplanet will have a compact envelope under certain con- 
ditions. What we cannot do with this model is obtain a 
quantitatively correct picture of such a compact envelope. 

Additionally, Fig. shows why the choice of EOS 
is not critical for the qualitative picture: Although the 
non-ideal effects might change the density stratifications 
of the compact inner parts, each solution which is not 
self-gravitating in its inner (barometric-law like) part, 
will have a counterpart solution with an effective core. 
Properties of the effective core will be dictated by the 
EOS, but its effect on the scale-height will remain the 
same. 

3.7. Estimating the Applicability of the Isothermal 
Assumption 

In the previous section we showed that an ideal gas is a 
good approximation for most of the parameter range we 
use. The validity of the isothermal assumption is examined 
below. 

By analogy with the pressure scale-height, a tempera- 
ture scale-height of a radiative stratification can be defined 
as: 



Hp 



'rad 



dr 
ainT 



(13) 
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Fig. 8. Cut through the envelope mass manifold, for a 
lO-^M^ core, a = 5.2 AU, and T = 123 K. Black circles 
represent the ideal gas, and squares are for the Carnahan- 
Starling EOS. 



where 



H 



dr 



P 



d\nP Q GM{r) 
for an ideal gas and hydrostatic equilibrium, and 



d\nT 



'rad 



ainP 



ZkLP 



(14) 



(15) 



47racGM(r)r4 ' 

"where a is the radiation constant, k is the gas opacity 
taken to be 0.1 m^ kg~^, c is the speed of light, and L is 
the core luminosity due to the planetesimal accretion rate 
of 10-6 Meyr-^ 

The temperature scale-height corresponds to the length- 
scale of a radiative giant-protoplanet over "which the en- 
velope temperature drops by a factor of 1/e. The spe- 
cific temperature scale-height HT{r)/rm\\ evaluates the 
ratio of the thermal length-scale to the radial extent 
of the entire protoplanetary envelope, at a position r. 
Evaluated at r = rniii, i?T(^Hiu) is the global estimate 
of the thermal scale- height of the protoplanet. Figure |51 
sho"ws that the isothermal assumption is valid for large 
portions of the manifold regions III and IV (cf. Fig. \^ , 
where HT(rHiii)/i'Hiii has values much larger than unity. 
These envelopes have a specific thermal scale-height above 
unity for at least the outer 90% of the envelope. Therefore, 
even though the small innermost envelope region is prob- 
ably non-isothermal, the protoplanet should be well rep- 
resented by the isothermal gas. 

Close to the giant-protoplanet 's critical core mass (e.g. 
logMcorc = —1.25 in Fig.0J, Hx/rHiu is expected to be of 
order unity and the isothermal assumption breaks down. 
Compact solutions (regions I and II from Fig. |21 and 
high Qcsg solutions in Fig.|5J) have very large Hx/rHiib in- 
dicating that nearly-vacuum space around the compact 
envelope is nearly isothermal. Detachment from a pro- 
toplanetary nebula could represent either hydrodynami- 
cally active protoplanets, or the collapsed gas giants with 



cleared protoplanetary nebula (i.e. mature giant planets). 
In both cases objects are expected to be deep in the non- 
isothermal regime. The radial profiles of the compact ob- 
jects will change if a detailed energy transport is included, 
but they will nevertheless remain compact. A compari- 
son of Jupiter's radius with that of our model planet (of 
equivalent mass and Tenv of 5000 K, estimated to be rep- 
resentative of Jupiter's average temperature from Guillot 
119991 shows that, with rcompact = 6.63 10^ m, our model 
falls short by less than 10% of reproducing the radius of 
the real gas giant. 

In the context of Jupiter's potentially rapid formation (or- 
der of 10^ years), it could be argued that the core accre- 
tion rate should be even higher. However, Hx/rHiU is pro- 
portional to the inverse of Mcore, and even if it is set to 
lO-'"^ M0yr-^, the validity of the isothermal assumption 
is still appropriate for the regions III and IV of Fig. |21 
Indeed, such high core accretion rates are applicable for 
cores comparable to (i.e. cores at late stages of a giant 
protoplanet's evolution), and are surely an overestimation 
for the younger cores (e.g. for the cores of 10"'^ Mq), mak- 
ing the case for the isothermal regime even more solid. 
However, because of the simplicity of our model, the re- 
sults are only qualitative, while quantitatively correct val- 
ues would only be accessible through a more elaborate 
model. 

Hx/rHiu shows that close to the critical core mass there 
are non-isothermal effects. 

But the basic isothermal picture is valid for most of (the 
quasi-homogenous part of) the manifold. It even appears 
that the possible transition from homogenous to compact 
state can be initiated within the isothermal regime. 



3.8. Manifolds and Environment 



Manifold solutions are dependent on four environmental 
parameters; the gas temperature T of the protoplanet (and 
of the surrounding- nebula), the orbital distance a from 
the parent star, the mean molecular weight /i, and the 
mass of the parent star Mi,. These parameters infiuence 
the balance of the two forces that determine the radial 
density structure - the outward force arising from the gas 
pressure, and the inward force of gravity; T and /x are 
connected with pressure through Eq. (^J, while a and 
determine the Hill-sphere, i.e. the volume of the envelope 
mass. 

Because of the simplicity of the model, the impacts of T 
and ^ on the solutions will be discussed together, as will 
the influence of a and M^. In reality, these parameters will 
have very different impact. 

Unless otherwise specified, the reference parameters 
throughout the current section are: log(Mcoro/^^e)="5, 
a = 5.2 AU, Tenv = 123 K, and ^ = 2.3 lO^^ kg mor^. 
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Fig. 9. Specific temperature scale-height as a function of 
the density at the core surface, for different subcritical core 
masses. Protoplanetary models with cores of -8 (black cir- 
cles), -5 (stars), and -3 (crosses) in logarithmic Mq units 
have HT(rHiii)/rHiii much larger than unity. This justifies 
the isothermal assumption for the manifold regions III and 
IV. 
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Fig. 11. Envelope mass solutions as a function of gas den- 
sity at the core surface, for orbital distances of 0.05, 0.1, 1, 
5.2, and 30 AU. Enlargement: The transition from uniform 
to compact envelope solutions is more abrupt for proto- 
planets at large orbital radii. This is a consequence of the 
larger Hill-sphere of outer protoplanets. 
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Fig. 10. Envelope mass solutions as a function of gas den- 
sity at the core surface, for gas temperatures of 100, 500, 
1000, 5000, and 10000 K. A change of T has no influence 
on the envelope mass of the non-self-gravitating regions, 
■while the same change of T -will produce a significant effect 
for protoplanets in self-gravitating regions. 



3.8.1. Temperature and Mean Molecular Weight 

Although this model is isothermal, the choice of gas tem- 
perature influences the solution manifold quantitatively. 
From Eq. Q it is clear that pressure relates linearly to 
temperature. Since the pressure force counterbalances the 



gravitational force, protoplanets "with hotter envelopes re- 
quire more gravity (and thus more mass) for a hydrostatic 
solution. The value of the critical core mass is a good ex- 
ample of the quantitative influence of the temperature. For 
example, the critical core mass for a 123 K protoplanet in 
Jupiter's orbit is 0.0948 M^, "while the critical core mass 
value for a 5000 K case is 24.5 M®. 

Figure E| sho"ws that, for subcritical cores and lo-w gas 
densities at the core surface (region III in Fig.|21l, the gas 
temperature has virtually no impact on the envelope mass. 
Since the envelope mass is small compared to the core 
mass, the envelope parameters (e.g. Tonv) have no influ- 
ence on the hydrostatic force balance via gravity feedback. 
On the contrary, for envelopes in "which self-gravity shapes 
the radial structure (regions IV and II in Fig.|21), the en- 
velope mass is significantly affected by different Tonv 
The scaling la"w "which relates manifolds of various tem- 
peratures is discussed in Sect. I3.1?1l 

As previously mentioned, this simple model does not in- 
corporate an energy transport equation, nor does it take 
into account the gas and dust opacities. Therefore, a 
change in /i cannot be distinguished from the correspond- 
ing change in T, and "will not be further discussed. 

3.8.2. Orbital Distance and Star Class 

Orbital distance, together -with the masses of the proto- 
planet and the parent star determine the protoplanet's 
gravitational sphere of influence, the so called Hill-sphere. 
Since the available volume for the protoplanet's envelope 
scales "with the cube of the orbital distance (see Eq. Q|), 
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the strength of the envelope's self-gravitating effect de- 
pends critically on the distance from the core to the parent 
star (see Fig. lllf) . Therefore, for the inner protoplanets to 
have (at least partly) self-gravitating envelopes, the gas 
density at the core surface must be larger than for the 
corresponding outer protoplanets. 

For solutions with compact envelopes (right side of Fig.lTTl 
and enlargement) the orbital distance has no impact on 
the envelope mass, since the radii of the compact inner 
part are typically several orders of magnitude smaller than 
their respective Hill-spheres. 

The transition from a uniform self-gravitating to a com- 
pact envelope is characterized by a considerable drop in 
the gas density for the outer envelope stratifications. In 
addition, protoplanets close to the parent star have rela- 
tively small Hill-radii and most of the envelope mass can 
be found in the proximity of the core. Therefore, the tran- 
sition from uniform to compact envelope for protoplanets 
close to the parent star is less abrupt than for more dis- 
tant protoplanets, as can be seen in the enlargement of 
Fig.lIT] 

Varying the mass of the parent star is equivalent to chang- 
ing the orbital distance of the protoplanet, provided that 
the gas temperature stays the same. It follows from Eq. 
that 6a~^ = SMi,, e.g. changing the orbital distance of the 
protoplanet from 5.2 AU to 1 AU is equivalent to chang- 
ing the mass of the parent star from iW^ = 0.21 Mq to 
Mi, = 30 M0. It remains to be seen whether this equiv- 
alence will hold for a more complex model, because the 
nebula properties will likely change in accordance with 
the known mass-luminosity relation as is varied. 

3.9. Static Critical Core Mass 

There are several definitions of the critical core mass cur- 
rently in use. The critical core mass concept has been in- 
troduced by various investigators (e.g. Perri & Cameron 
('1974'), Mizuno et al. l(TU75jl . Mizuno H1980|l . Bodenheimer 
& Pollack l(TU5Hjl . Wuchterl H1991a|l ). As a starting point, 
we choose here a definition suggested by Wuchterl H1991afl . 
for 'static critical core mass': No more static core- envelope 
models with increasing core mass exist at the critical mass. 
This definition is valid along a (time) sequence of proto- 
planetary models with increasing Afcore- It is only along 
such a sequence, in the context of the static models, 
that a time evolution with growing cores can proceed. 
Essentially, the static critical core mass is the largest core 
mass for a static protoplanet that can be embedded in 
a given nebula, characterized by a nebula gas density, a 
temperature, and a distance from a parent star. 
For the (a = 5.2 AU and T = 123 K) manifold 
this means that, among the solutions with £ionv(''Hiu) — 
1.4 10~® kgm~^ (defined for the minimum mass solar neb- 
ula, e.g. Hayashi et al. ll985|l . the solution with the largest 
core mass determines the static critical mass (Fig.^J the 
innermost solid line). This gives a static critical core mass 
ofM.^.'^SN = 0.0948Me. 



Figure IT^ shows that the value for the critical core mass 
exhibits a generally weak dependence on the density of the 
surrounding nebula, so the choice of from different 

nebula models is not critical. For the very dense nebulae 
(around 10~^ kgm"'^) and depending on the choice of the 
solution branch (cf. Sect. I3.9"T|| . the values for the local 
critical core masses can span several orders of magnitude 
even for the same nebula. 

The critical core masses for different manifolds are pre- 
sented in Table El and are found to depend on the param- 
eters that affect the hydrostatic balance (cf. Sect. 13. 8|) . 

By comparing Figs. ^ and [51 it follows that the natural 
choice for the global static critical core mass, one which 
is valid for the whole manifold, should be the core of the 
protoplanet which is at the interface of all four manifold 
regions (cf. Fig. El). The model at the interface has a min- 
imum in the envelope mass, for a manifold cut along the 
constant Qcsg value. The interface is also an inflection point 
for a manifold cut at a constant Mcore- The conditions for 
the global static critical core mass thus are: 

(16) 



Since the numerical values for the global critical core 
masses are very close to the values of the critical core 
masses from the definition suggested by Wuchterl H1991a|l . 
we do not present the global numerical values separately. 

The values obtained for critical core masses in this 
model agree well with those of Sasaki H1989|l . who used a 
similar set of assumptions. However, such isothermal val- 
ues are significantly smaller than today's commonly ac- 
cepted critical mass values, obtained with the inclusion 
of detailed energy transfer, which are typically between 7 
and 15 Mg. The reasons for this are two- fold. Firstly, we 
use the equation of state for an ideal gas. Secondly, the 
temperature of the isothermal gas is taken from nebula 
models, hence the nebula temperature is the temperature 
of the entire protoplanet. This is certainly a lower limit for 
a realistic temperature of the interior of the protoplanet. 
Larger critical core mass values are obtained if the gas 
temperature is in the range of the temperatures for the 
interior of gas giants modelled with detailed energy trans- 
fer (cf. Sect. 13.8.1)) . Clearly, the correct determination of 
the critical core mass requires temperature structure, but 
the emphasis in this work was not on quantitative details, 
but rather on global qualitative features. 

3.9.1. Local Critical Core Mass 

From Figs . 1 1 21 and 1131 one can see that, for each subcritical 
core immersed in a nebula, there are at least two permitted 
solutions. However, if one considers only the time-sequence 
of hydrostatic models with a growing core, it is clear that 
solutions with higher density at the core surface cannot 
be reached. 
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Fig. 12. Solution branches - isobars for gonv(''Hiii) — Pout 
- for (a = 5.2 AU and T = 123 K) manifold: the stan- 
dard solar nebula solution branch is represented by the 
innermost solid line; an enhanced nebula with gout — 
10~^kgm^'^ (dashed lines) has multiple solution branches; 
each solution branch has its own maximum core mass, and 
hence local critical mass 



The situation is more complicated if the protoplanet is 
embedded in a denser protoplanetary nebula. Our model 
clearly predicts multiple solutions for certain sets of pa- 
rameters (Fig. 113(1 . Instead of one solution branch for a 
given nebula (with two solutions for each core, as for a 
minimum mass solar nebula), several solution- branches 
are possible, again each with two solutions for a spe- 
cific core (Fig. 1121 dashed solution branches for gout = 
10~^ kgm""^). Multiple solution-branches are enabled by 
envelope self-gravity (cf. Fig. I14|l and are due to tidal re- 
strictions imposed by the parent star via rniii (cf. Fig. O 
region IV). 

Each solution branch has one critical core mass, be- 
yond which there is no static solution, for a sequence of 
hydrostatic models with increasing core mass. For the min- 
imum mass solar nebula this means one critical core mass, 
according to the definition suggested by Wuchterl ()1991a|l . 
For some denser nebulae, however, the existence of several 
branches implies several local critical core masses, where 
solutions beyond the critical core mass of the branch are 
unavailable locally. After reaching the local critical core 
mass, the planet could, in principle, continue evolving by 
'jumping' to another branch. One of us has observed sim- 
ilar behavior for certain sets of initial parameters in hy- 
drodynamical models. The local critical core mass satisfies 
the above definition but not Eq. ((16|l for the global critical 
core mass. 
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Fig. 13. For nebula density enhanced relative to a 
minimum-mass solar nebula, even more than two hy- 
drostatic equilibria could exist; M: protoplanetary solu- 
tions with logMcorc/[AfEarth] = -2 that fit into gout = 
lO^^kgm^^ nebula; DS: double solutions, a special case 
of multiple solutions, cf. Figs. 131 and 1 151 S: protoplanetary 
solutions with the same core, whose envelope fits into the 
minimum-mass solar nebula. 
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Fig. 14. Density profiles for the solutions which fit into the 
same (10^® kgm^'^) nebula. These solutions are labelled 
with M in Fig. ^ 



3.9.2. Double Maxima 

A special case of multiple solutions can be seen in Figs.|31 
1131 and 1 151 as double peaks in the envelope mass. For ev- 
ery (subcritical) core, two special solutions, which fit into 
the same nebula cloud (i.e. have the same f?(rHiii)) and 
have almost exactly the same envelope mass (equal to 
one part in lO'*, or better), are found to exist. Usually 
these two solutions have a very similar stratification in 
the outer parts of the protoplanet's envelope, but deep in- 
side the protoplanet their radial structure is quite different 
(cf. Sect.ES). 
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Fig. 15. Mass and density radial structure of the special 
case of multiple solutions, where two protoplanets have 
the same core, almost the same envelope mass, connect 
to the same nebula, but have different radial structure. 
These solutions are labelled DS in Fig. El 



Supercritical cores do not posses such a feature, be- 
cause the density profile always effectively goes to zero 
long before the Hill radius is reached. Therefore there is no 
significant contribution to the envelope mass in the outer 
stratifications, and the envelope mass increases monoton- 
ically with the gas density at the core surface (cf. Fig.|2I). 

Envelopes with lower gas density at the core surface, 
£)J.°g, (Fig. 115(1 have a maximum possible mass (for the 
corresponding manifold) because the envelope gas density 
at the core surface is low enough to ensure uniformity for 
the major part of the radial structure. Consequently, the 
envelope density does not substantially decrease from the 
core-surface value. At the same time, Qcsg is high enough 
to allow significant mass contributions from the outer 
parts of the envelope, where the volume (and therefore the 
mass, for a given density) per unit radius, is the largest. 
Values for such maximum envelope masses are tabulated 
in Table El 

Envelopes with higher g^^^ build up the self-gravitating 
effect (SG effect starts as soon as Monv ~ ^^corc) very close 
to their core. Because of the very strong self-gravitating 
effect (Monv ~ SMcoro for the innermost regions), the ra- 
dial density fall-off close to the core (Fig. I15|l is strong. 
A new, effective core is formed from the dense envelope- 
gas-layer wrapped tightly around the core. In this case, 
the envelope density distribution resembles one with the 
core (and the radius) of this effective core. In stratifica- 
tions where the envelope mass becomes comparable to the 
effective core, another self-gravitating effect changes the 
radial envelope density distribution to gonv oc r~^. 
For a particular choice of £ics|'^, the envelope density pro- 
file in the outer stratifications matches that of thus 
making the mass of both envelopes almost equal. 



3.10. Temperature-Mass Invariance 

It has been noted that, if mass and distance are mea- 
sured in a system of appropriate units (i.e. mass in units 
of core mass, and distance in units of core radii), solution 
manifolds with different temperatures are almost identi- 
cal, except for a shifting on a core-mass-axis, according to 
the relation: 



1l = (Ml)'^/3 

T2 ^Nh' 



(17) 



that can be derived for homologous envelopes satisfying 
£'i(''i/?'coro,i) = Q2[r2/rco'cc,2), for any pair of ri and r2 
such that ri/rcorc.i = f2/fcorc.2- In other words, the ra- 
dial profile of a certain protoplanet with core mass Mi 
and temperature Ti will be the same as the radial profile 
of another protoplanet with core mass M2 and tempera- 
ture T2, if Eq. 1(17(1 is obeyed, and if the mass is measured 
in units of core mass and the length in units of core radius. 
This is true for all manifold regions, sub- and super- 
critical, self-gravitating or not. Note that in Fig. 1101 the 
non-self-gravitating region was not affected by a change in 
envelope temperature, but relation l(17() does hold even for 
non-self-gravitating envelopes, since it connects envelopes 
with different temperatures and core masses. Fig. llUl was 
plotted for different temperatures, but constant core mass. 

4. Discussion and Conclusions 

In an effort to obtain a global overview of hydrostatic pro- 
toplanetary equilibria, we have chosen a simple physical 
model so as to be able to clearly understand the interac- 
tion of competing processes. 

Our use of relatively simple physics has several conse- 
quences; because the ideal gas equation of state is used, 
gas particles are 'soft', and can be compressed as much 
as is needed, in effect overestimating the importance of 
gravity relative to gas pressure, when large envelope-gas- 
pressure is applied. A comparison of the ideal gas EOS 
to the numerical Saumon-Chabricr EOS shows disagree- 
ment for the log T = 2.1 isotherm and densities above 
Q = AO kgm^"^. This would indicate that the non- ideal 
EOS is needed for high-density effective cores. 

It has been noted that manifold properties are insensi- 
tive to variation of orbital distance a or mass of the parent 
star Mi,, as long as aM~^ = const holds (cf. Sect. 13. 8T^ . 
Also, solutions whose envelope temperature and core mass 
obey relation ((17(1 are found to be the same, if appropriate 
units for mass (i.e. core mass) and length (i.e. core radius) 
are used. This indicates the existence of analytic solutions 
for some envelope regimes, through certain dimensionless 
scaling variables. Such a treatment is, however, out of the 
scope of the present paper. 

An envelope gas temperature is equal to the nebula T 
throughout the protoplanet, and that certainly underes- 
timates the thermal pressure and hence reduces the val- 
ues for the critical core mass. However, from Eau. ll7( one 
can show that for a more realistic estimate of the envelope 
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temperature representative for the young planets (5000 K) 
critical core mass values are overestimated (~ 24. SM^), 
because of envelope isothermality/lack of an energy trans- 
port equation and use of ideal-gas EOS, when compared 
to canonical critical core mass values from protoplanetary 
models with detailed microphysics. 

Both the local and the global critical core masses signal 
the end of the availability of the hydrostatic solutions. In 
the case of the local critical core mass, non-availability 
holds for a small region of the parameter space around 
the local critical core mass, while for the global critical 
core mass this is true for every core larger than the critical 
core mass. The significant difference between the two types 
of critical core mass is that, at the global critical core 
mass (and above), the non-isothermal effects are crucial 
in shaping the structure of the protoplanetary envelopes, 
and are present throughout the parameter space. These 
non-isothermal effects are important for determining the 
details of the dynamical disk-planet interaction. 

The critical core mass values obtained in this model are 
almost two orders of magnitude smaller than the canoni- 
cal critical core masses which incorporate detailed energy 
transfer. Thus, if subcritical or just-critical regimes of a 
dynamical disk-planet interaction are to be investigated 
through a model that is locally isothermal, the planet mass 
should be set appropriately. Most of the present locally- 
isothermal disk-planet models (e.g. Kley 1999, D'Angclo 
et al.EnnSEnni Nelson & PaDaloizou l2004,) operate with 
planets which should be deep in the super- critical regime. 

A solution set from our model encompasses solutions 
that are reminiscent of the planets in the various stages of 
evolution (from small rocks embedded in the dilute nebula 
to the mature planets as we know them) , and of various 
configurations (the telluric planets of region I in Fig. |21 
and the gas giants of region II). The 'nebula' and 'ma- 
ture planet' regimes are the physically intuitive beginning 
and end phase of planetary evolution. However, the 'pro- 
toplanet' regime presents us with an interesting region in 
parameter space, where planets could make the transition 
from 'infancy' to 'maturity'. Depending on the detailed 
structure and the dynamics of the surrounding nebula, it 
is easy to conceive a standard scenario of planet formation. 
That is, the accretion of nebula gas onto a supercritical 
protoplanet. Other scenarios could be imagined as well, 
e.g. a massive protoplanet could release a major part of 
its envelope to reach the appropriate equilibrium, or it 
could dramatically condense its otherwise mostly gaseous 
envelope. Amounts of dust in the environment will doubt- 
less play a very important role in the process. 

In conclusion, several important features of the solu- 
tion set have to be mentioned: 

1. Two basic types of the envelope equilibria are found 
for protoplanets: 

• uniform; the density of the envelope gas drops 
weakly from the core to the outer boundary 

• compact; a dense gas layer forming an effective 
core, and a very low, exponentially decreasing, gas den- 



sity further out 

Both types can be self-gravitating or non-self- 
gravitating, dividing the solution manifold into four 
distinct regions. 

2. As a consequence of the envelope's self-gravitating ef- 
fect, a wide range of possible envelope solutions exists. 

3. We have developed a new concept for the global static 
critical core mass, which marks the contact point of 
all four qualitatively different types of protoplanets. 
This concept is based on a qualitative change of the 
envelope properties while considering a complete set of 
available solutions (a solution manifold), as opposed to 
the critical core mass definitions which are valid only 
for a solution subset fitting a particular nebula. 

4. For every subcritical core there are at least two en- 
velope solutions possible (a self-gravitating one and a 
non-self-gravitating one) for a given nebula, and for a 
certain nebula parameters the number of the possible 
envelope solutions can be even larger. Such nebulae 
also have multiple (local) critical core masses. 

5. The global static critical core mass value is shown to 
decrease with the increasing orbital distance a, mainly 
because of the decrease in the temperature of the sur- 
rounding nebula. 
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